Phase velocity and particle injection 
in a self-modulated proton-driven plasma wakefield accelerator 
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It is demonstrated that the performance of the self-modulated proton driver plasma wakefield 
accelerator (SM-PDPWA) is strongly affected by the reduced phase velocity of the plasma wave. 
Using analytical theory and particle-in-cell simulations, we show that the reduction is largest during 
the linear stage of self-modulation. As the instability nonlinearly saturates, the phase velocity 
approaches that of the driver. The deleterious effects of the wake's dynamics on the maximum energy 
gain of accelerated electrons can be avoided using side-injections of electrons, or by controlling the 
wake's phase velocity by smooth plasma density gradients. 
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A plasma is a promising medium for high gradient ac- 
celeration of charged particles. It can sustain fields orders 
of magnitude higher than the breakdown fields of con- 
ventional accelerators |T]. One can excite strong plasma 
waves either with lasers or with charged particle beams 
[H E]. One of the very attractive approaches is to use 
already existing TeV proton beams as a driver to gener- 
ate plasma wake fields. Due to the limitation set by the 
transformer ratio, the energy gain of the witness beam 
cannot be much larger than the driver energy Em- 
ploying a TeV proton driver allows in principle to acceler- 
ate an electron bunch to TeV energies in one single stage 
thus alleviating the technical burden of multistaging. 

It has recently been shown using detailed simulations 
O E] that a high gradient plasma wake fields can be 
generated with an ultra-short bunch of protons. In 
that scenario, the proton bunch was shorter than the 
plasma wavelength. Unfortunately, such ultra-short pro- 
ton bunches are not presently available. The length of ex- 
isting TeV-class proton bunches is of order Laz ~ 10cm. 
The characteristic plasma field, the so called wave break- 
ing field is ii'wB = mcwp/e « 96y^7i7(cni~^V/m, where 
(jjp — \J Airnee^ / Ta is the plasma frequency defined by the 
electron density Ue. Accelerating gradients of a GV/m- 
scale require a plasma density of at least rig — lO^'^cm^'^ 
corresponding to the plasma wavelength Ap = iir/kp = 
2TTc/ujp ^ 3nim. Thus, the existing proton bunches cor- 
respond to (Tz = (10 — 100)Ap and cannot efficiently gen- 
erate wake fields in such plasma. The situation with 
the proton bunches is very much the same as it was 
with laser pulses in the 1980's. The availability of long 
laser pulses necessitated the invention of a self-modulated 
laser wakefield accelerator (SM-LWFA) |7] . Subsequent 
progress in ultra-short pulse laser technology removed 
the need for self-modulation and led to successful mono- 
energetic electron acceleration in the bubble regime [8j 
that reached GeV energies. 



A long proton bunch propagating in an overdense 
plasma is also subject to self-modulation at the back- 
ground plasma wavelength [9\. The effect of self- 
modulation opens a possibility to use existing proton 
bunches for large amplitude wake field excitation. An 
experimental program is currently under consideration 
at CERN. The Super Proton Synchrotron (SPS) bunch 
with 450 GeV protons is proposed as the driver for the ini- 
tial stage of the experimental campaign. The wake field 
will be used for accelerating externally injected electrons. 
The injected particles must be trapped in the wake field. 
The trapping condition depends on the wake field ampli- 
tude and phase velocity [1^ . Because it is expected that 
the SPS bunch will generate a weakly-nonlinear plasma 
wave with the same phase velocity Wph as the speed of the 
driver, it is natural to assume that the gamma-factor of 
the injected electrons 7ci must be comparable to that 
of the proton driver 7p for them to be trapped. As 
demonstrated below, that is not the case because the 
spatio-temporal nature of the self-focusing instability of 
the proton bunch considerably reduces Vp^. 

Although it has been realized for some time [TT] that 
the phase velocity of the plasma wake produced by the 
self-modulation instability of a laser pulse is slower than 
the pulse's group velocity, this was not an important is- 
sue because the laser group velocity was usually modest. 
For the self-modulated proton-driven plasma wakefield 
accelerator (SM-PDPWA), the wake slowdown is of crit- 
ical importance. Here we show that the phase velocity of 
the unstable wave is defined not so much by the driver 
velocity, but mainly by the instability growth rate. The 
wake field is greatly slowed down at the linear instability 
stage when the growth rate is at its maximum. At the 
nonlinear saturation stage, the wake reaches the driver 
phase velocity. We also propose a method to manipulate 
the wake phase velocity by smooth longitudinal density 
gradients. 
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To describe the wake slowdown analytically, we adopt 
the formalism developed within the framework of the 
envelope description of the driver j9] . We assume an ax- 
isymmetric bunch driver and utilize the co-moving and 
propagation distance variables ^ = f3oct — z and t — z/c, 
respectively, where /3o = Vz/c (vz is the velocity of the 
bunch) and z is the bunch propagation direction is z. 
Further, the driver bunch is assumed stiff enough so that 
its evolution time is slow dr <C cd^. The bunch is as- 
sumed to be long: ujpaz/c 3> 1. The assumed bunch 
density profile is p(r, ^) = pQ ipir)f{^), where po = rihe is 
the charge density of the proton bunch. For simplicity, 
the step-like radial profile ■(/'(r) = 6(rb — r) is assumed, 
where rb(^, r) is the evolving radius of the bunch's enve- 
lope, and 0(rf, — r) is the Heaviside step-function. The 
betatron frequency and wavenumber of the self-focused 
bunch are defined as copo = ckpQ ~ y^Annbe'^ /2^bmp, 
where nip is the proton mass. In the limit of a thin bunch 
(kpTb 1) and linear plasma response, the equation of 
motion for the bunch's radius (in normalized coordinates 
h = n/no, T = i^iioT: C = fcpC) is given by [5] 
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^sin(^-e')de', (1) 
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where in = ere/ (kporbo) with ae, rbo being the bunch an- 
gular divergence, initial radius, and longitudinal current 
profile, respectively, and !{£,) = fiO^bi^'^)- Perturbing 
Eq.Q about the initial radius f ^ = 1 + ^^b, yields the 
linearized equation 



d^Srb 



+ 3e„5rb = Sn. 
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Following the approach of Bers [12_, we find an asymp- 
totic solution of this equation for sufficiently late times, 
f > Lg, where ^ l/F is the e-folding length, and 
F is the growth rate of the instability. Substituting of 
drb = Sfb exp{—i6ujT + iki) into Eq. (jij) yields the dis- 
persion equation 6ijj{k). The peak growth rate is calcu- 
lated [12' by introducing — Su ~ vk, where v = ^/f , 
and requiring that D{duj ,k) — and dD{Suj , k) /dk = 0: 

(l-fc2)M..2 + 3el) = 1, (3) 
-k{~Suj^ +3il)- 6ujv{l- k^) = 0. (4) 

Equations (|3| and Q lead to the standard dispersion 
relation typical of the beam-breakup instability which 
is known [17 \ to always possess a growing mode with 
> 0. To simplify the algebra, we assume that 3e^j 
is small and consider the initial stage of the instability 
corresponding to fc = 1 + 6k, where 6k <ti 1. In this 
limit 6uj ^ €n, and Eqs. ([3|,(|4]) reduce to the dispersion 
relations 6uj^ = —v, and 6uj — ~2v6k. The complex roots 
are given by 



6uj+ = 



:\ 1/3 



In physical units, the condition for 6uj ^ e„ can be ex- 
pressed as 



_ 2^ ribm 

t ^ ishort = 7fc 
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where e„ = "/bCgrb is the normalized beam emittance. 
The number of e-foldings is given by = '^{6uj)f — 
^{6k)£^ = 33((5cL>)f/2, and therefore the growth rate 

F = 3'i{6uj) 12 = 3^3 (|/(|) /f ) ^ ^ /4. In dimensional 
variables, the instability growth rate is expressed as 
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and the maximum number of e-foldings achieved at t — 
ishort is iVe « (3y3/4)|i/3f2/3 ^ {3V3/A)kpazin^^^. 
Note that for the typical parameters of the SPS bunch 
and plasma density Ue = 10^* . . . lO-'^^cm"'', ~ 100. 
Therefore, it is most likely that the self-modulation in- 
stability will enter the nonlinear regime prior to t = tghort , 
and the above assumptions will remain valid throughout 
the linear stage of the instability. 

The crucial observation is that 6lu and 6k have not 
only imaginary parts responsible for instability growth, 
but also real parts. It is these real parts that change 
the wake phase velocity. The wake phase is ipe = (1 + 
^{6k))£^ — ^{5uj)t, and the phase velocity Uph = Vb + 
{dip/df)/{d'ilj/dS,)ujpo/u}p. Substituting (|5| for 6uj and 
6k and neglecting small terms on the order of v^^^^, we 
obtain the phase velocity of the growing mode 
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The wake phase velocity (|8| can be significantly lower 
than the speed of the bunch due to the instability dis- 
persion. The relativistic 7— factor of the wake phase ve- 

/ N-i/2 
locity 7pft = ( 1 — v^y^jc? \ can be an order of mag- 
nitude lower than that of the driving bunch. This effect 
will prohibit electron acceleration to high energies at the 
growing instability stage. Yet, one can easily see from 
the formulas that the phase velocity decrease is closely 
connected to the instability growth rate. Thus, one may 
expect that when the instability saturates, the phase ve- 
locity of the wake becomes close to that of the bunch. 
This effect might help to inject low energy electrons into 
the wake of a highly relativistic proton bunch at a later 
stage of the instability, just before the nonlinear satura- 
tion of the instability. 

We should mention here that the dispersion relation 
allows also for a purely oscillating mode with 6loq — 

/. \l/3 / .\2/3 

— (C/'t) , 6k = \ ■ This mode has the su- 
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Figure 1: (color online) A 3D hybrid PIC simulation of 
SPS bunch self-modulation in plasma with constant density 
rie = 7.76 ■ 10^* cm~^. Frame (a) shows the evolution of the 
maximum accelerating field. The linear instability persists 
for some 5 meters before the nonlinear saturation. Frame (b) 
shows the phase velocity of the wake as a surface dependent 
on the propagation distance z and the bunch coordinate ^. 
The phase velocity experiences a deep drop towards the tail 
of the bunch at the linear stage of the instability. When the 
bunch is completely modulated and the instability saturates 
at z > 5m, the wake phase velocity surges up to that of the 
driving bunch. Frame (c) gives a snapshot of the wake phase 
velocity after 2.5 meters of propagation distance (the solid red 
line). The broken line gives the analytical expression Q. 



perluminous phase velocity 
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Stable propagation of a modulated beam has been ob- 
served recently in simulations |13] . 

The envisioned experimental program at CERN will 
use the SPS bunch. It normally delivers 1.15-10^^ protons 
at 450 GeV/c with the normalized emittance e„ = S.S/um 
and the length = 12cm. We use these bunch parame- 
ters in our 3D PIC simulations with the newly developed 
hybrid code H-VLPL3D [T3]. This new code simulates 
the background plasma hydrodynamically while high en- 
ergy bunches are treated with a full kinetic algorithm. 
The hydrodynamic part of the code introduces much less 
numerical dispersion into the plasma waves than a PIC 
code with the same resolution. The bunch focused to 
CTr = 0.19 mm is sent through plasma with the free elec- 
tron density rig = 7.76- lO^^'cm"'^. The maximum bunch 
density on axis is rif, = 1.5 • lO^^cm"^. To avoid beam 
hosing and to seed the self-modulation, we assumed the 
bunch is hard-cut in the middle [S] . 

The simulation results are presented in Fig. [T] From 
the frame (a) we see that the linear instability stage 
lasts during the first 5 meters of propagation. Then, the 
bunch is completely modulated and the nonlinear satu- 
ration is reached. The wake phase velocity is shown in 
the frame (b) . A significant slowdown of the wake is ob- 
served during the instability and along the bunch. The 
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Figure 2: (color online) Electron energy spectra nt z = 10 m. 
There were two groups of electrons. One was injected on- 
axis at the plasma entrance; another was side-injected at the 
angle of 0.005 radian. In the latter case, the electron bunch 
trajectory should cross the driver axis at z = 6m, ^ = 20cm 
behind the bunch head. The on-axis injection lead to a broad 
spectrum and low energy gain. The side injection resulted in 
a quasi-monoenergetic beam at 1.2 GeV energy and about 1% 
energy spread. Initial electron energy was 20 MeV for on-axis 
injection and 10 MeV for side injection. 



lowest phase velocity is The lowest relative phase veloc- 
ity is (uph — c)/c w —3.5 • 10^^, corresponding to a wake 
Iph ~ 38. This is an order of magnitude lower than the 
7— factor of the driving bunch. Frame (c) compares the 
simulation result (the solid red curve) with the analytic 
expression ([8|, in which we substituted the SPS bunch 
parameters. This snapshot of the wake phase velocity 
has been taken at z ^ 2.5 m, in the middle of the linear 
instability stage. A reasonable agreement between the 
simulation and the analytical theory is observed. 

The wake slowdown has a dramatic impact on the elec- 
tron trapping and acceleration. First, it allows for trap- 
ping of low energy electrons whose velocities are compa- 
rable with the wake velocity. However, the energy gain in 
the slow wake is very limited due to fast dephasing. The 
energy gain is given by AW « 'ij^h'^c^i^max/EwB) 
|10| . At the linear stage we have Emax ^ E\yb and 
the energy gain is low for small "fpfi. The dephasing, 
however, has a much worse effect : if the electrons have 
been injected into the early instability phase of the slow 
wake, they can be lost when they overtake the wave and 
enter its defocusing phase. The dephasing distance is 



27r7p^ and for the slow wake field it can be shorter 
than the distance needed for the instability to develop. 
For this reason, the electrons must be injected late in the 
development of the instability, when the phase velocity 
begins to grow. In our simulation, the optimum point 
for injection is located around z = 5m. The wake phase 
velocity is still low here, but starts growing rapidly as the 
bunch reaches complete modulation. 

A possibility to inject electrons into the wake is side 
injection [TS]. In this bunch of electrons is prop- 

agating at a small angle with respect to the driver. The 
advantage of side injection over on-axis injection is that 
electrons are gradually "sucked-in" at the right phase by 
the wake transverse field. This leads to high quality 
quasi-monoenergetic acceleration of electrons. 
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Figure 3: (color online) Test electron bunch (dark green discs) 
and the accelerating field (the red-blue wave) at 2 = 5 m. The 
electrons are captured and split into micro bunches located 
exactly in the accelerating and focusing phases of the wake. 



In our simulation, we have injected two bunches of test 
electrons. The first bunch of 20 MeV electrons was in- 
jected on-axis at the plasma entrance. This electron en- 
ergy roughly corresponded to the minimum wake phase 
velocity at the tail of the driver. We found that dur- 
ing the linearly growing instability stage, these electrons 
underwent more than one oscillation in the ponderomo- 
tive bucket. Finally, after 10 meters of propagation, the 
maximum energy gain was about 200 MeV with a rather 
broad energy spectrum as seen in Fig. |2] 

The second electron beam with 10 MeV energy was 
side-injected at an angle of 0.005 radian. The electron 
bunch trajectory was designed to cross the driver axis at 
z = 6m, <^ = 20cm behind the bunch head. Due to the 
small injection angle, however, the electrons were sucked 
in into the wake much earlier, at the position z k, 5m. 
The wake transverse fields have put most of the bunch 
electrons into the focusing and accelerating phase. The 
electron beam and field configuration just after the elec- 
trons entered the wake is shown in Fig. [3j The electron 
beam is split into micro-bunches located exactly in the 
accelerating and focusing phases of the wake. Due to this 
configuration, the side injected beam resulted in a maxi- 
mum energy gain of 1.2 GeV and a rather narrow energy 
spectrum. 

The low energy spread and efficient acceleration of the 
side injected electrons are also due to the fast rise in the 
wake phase velocity just after the injection position, as 
seen m Fig. [l]^b). The electrons gain energy while staying 
in the accelerating phase of the wake. 

We have seen above that the growing mode (|5]) has 
the low phase velocity ([s]). There is a way, however, 
to manage the phase velocity of the wake by employ- 
ing a gentle longitudinal plasma density gradient as it 
was discussed in [16j. To elucidate the effect, we have 
performed an additional simulation with the same beam 
parameters, but introducing a positive plasma density 
gradient: ne(z) = neo(l + -z/d) with the characteristic 
length d = 200m. The phase velocity obtained in this 
simulation is shown in Fig. |4j 

The phase velocity at the head of the beam takes a dive 
as defined by the growing mode dispersion. However, 
the positive plasma density gradient compensates for the 
mode dispersion and at the tail of the beam the wake 
phase velocity becomes equal to the speed of light and 
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Figure 4: (color online) Phase velocity control by positive 
plasma density gradient, a) Phase velocity at z=2.5 m. b) 
Full dynamics picture; the transparent plane marks the speed 
of light. 



even superluminous. 

In summary, we have shown that the self-modulational 
instability of a charged beam in plasma corresponds to 
a growing mode with a slow phase velocity. The wake 
velocity is much lower than that of the driver. The wake 
slowdown is due to the real part of the frequency of the 
unstable mode. Although this effect limits electron en- 
ergy gain at the stage of the linear instability growth, the 
low phase velocity can be harnessed to inject low energy 
electrons into the wake of a highly relativistic driver. We 
also have shown that the side injection of electrons at a 
small angle with respect to the driver axis may drasti- 
cally improve the quality of acceleration. The transverse 
field of the wake sucks in the injected electrons and au- 
tomatically puts them into the right acceleration phase. 
Finally, we show that the wake phase velocity can be 
controlled by a longitudinal plasma density gradient. 
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